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ABSTRACT 

The fractional Brownian motion with index a is introduced to construct the fractional 
excursion set model. A new mass function with single parameter a is derived within 
the formalism, of which the Press-Schechter mass function (PS) is a special case when 
a = 1/2. Although the new mass function is computed assuming spherical collapse, 
comparison with the Sheth- Tormen fitting function (ST) shows that the new mass 
function of a » 0.435 agrees with ST remarkably well in high mass regime, while 
predicts more small mass halos than the ST but less than the PS. The index a is the 
Hurst exponent, which exact value in context of structure formation is modulated by 
properties of the smoothing window function and the shape of power spectrum. It is 
conjectured that halo merging rate and merging history in the fractional set theory 
might be imprinted with the interplay between halos at small scales and their large 
scale environment. And the mass function in high mass regime can be a good tool to 
detect the non-Gaussianity of the initial density fluctuation. 
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1 INTRODUCTION 

Halo models are widely applied in the campaign of cosmo- 
logical parameter estimation in precision from cosmic large 
scale structures as well as in the expedition of understand- 
ing structure formation. The mass function is a fundamen- 
tal ingredient of halo models. The most famous analytical 
formula of mass function, the P ress-Schechter mass func - 
tion (hereafter PS) was derived bv lPress fc Schech ter ( 1974) 
based on the spherical collapse model. The PS function can 
be alternatively deriv ed with the rand om walk or the ex- 
cursion set formalism l|Bond et al.lll99lf) . By smoothing the 
linear density field on different scales with a sharp k-space 
filter, the density fluctuation within the characteristic scale 
could be regarded as a random walk against the variance 
of the smoothed field at this scale. Consequently the whole 
theory of random walk can be grafted to model the den- 
sity contrast field. The elegant theory provides a concise 
analytical framework to study various processes in cosmic 
structure formation, and is embr aced with grea t inte rests 
by the community. For instance, lLacev fc Cold l|l993f) ex- 
plicitly calculated the merger rate, halo formation time, 
and re levant properties of galaxy clusters; ISheth fc Tormenl 
( 2002) adopted the excursion set theor y with moving barrier 
to study ellipsoidal collapse of halos: IZhang fc Hul l|200fil) 
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solved the excursion set theory with moving barrier of arbi- 
trary shape and discussed the HII bubble size during re- 
ionization; an d voids phenomenon is e xplored within the 
framework bv lFurlanetto fc Piranl j200d) . 

The success of the random walk formalism in cosmol- 
ogy is prominent, but the primary product of the excursion 
set theory, the PS mass f unction, is a poor description to 
simulations at all epochs jReed et aUl2006ft . The common 
practice is to parameterize the PS function, and then fit the 
function to simulations to pin down free parameters (e.g . 
Sheth fc Tormenl 19991: 1 Jenkins et all200ltlReed et all200a 
Warren et alJl200d) . Many functional forms have been pro- 
posed by various a uthors to account for di fferent effects of el- 
lipsoi dal collapse (ISheth fc Tormenll2002l) , angular momen- 
tum (DclPopolo 2006a) and the index of power spectrum 
jReed et al.ll2006MBetancort-Riio fc Montero-Dortal fcOOd) 
claims that the "all-mass-at-center" problem shall be prop- 
erly formulat ed to obta in the correct mass function in high 
mass regime, iLed (120061) assumes there is a break in the hi- 
erarchical merging process and obtains much shallower mass 
function in low mass regime. 

In this report, we construct a fractional excursion set 
theory by replacing the conventional random walk with the 
fractional Brownian motion of index a. The standard ex- 
cursion set theory is simply a special case of the the new 
theory. The difference between the normal random walk 
and the fractional random walk lies in that the latter takes 
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the correlation between walking steps into account. A new 
mass function is derived with the fractional excursion set 
theory, which contains one parameter a in connection with 
the correlation of steps of the random walk. Although the 
new mass function is derived with the boundary condition 
of a single fixed absorbing barrier, i.e. in spherical collapse 
scenar io, it is in good agreeme nt with the Sheth-Tormen for- 
mula dSheth fc Tormenlll999|. hereafter ST) with a « 0.435 
in high mass regime, while has more small mass halos than 
ST and less small mass halos than PS. 

The layout of this paper is that at first we recite the 
excursion set theory briefly in Section 2, then in Section 3 
we introduce the fractional Brownian motion to develop the 
fractional excursion set theory and subsequently derive a 
new mass function, and the last section is of discussion. 



2 THE EXCURSION SET THEORY 

The initial density fluctuation 8 — p/p~ l<JClin early uni- 
verse is Gaussian and evolves linearly. If the density contrast 
in a region exceeds a critical value 8 C , the mass in that region 
will collapse and be virialise d in future to form a halo. 

As pointed out bv lBond et alJ dl99ll) . at an arbitrary 
point in the universe, the density contrast smoothed with a 
window function Wm(R) of characteristic scale R is a func- 
tion of the underlying total mass M(R) ~ pR 3 included by 
the smoothing window, the M effectively represents the scale 
R. The variation of the smoothed density contrast S(M) 
forms a trajectory in the plane of S(M)-M. The collapsing 
condition S c is turned into an absorbing barrier over the 
trajectory, at the largest M where S(M) firstly crosses the 
barrier, the trajectory will be absorbed, i.e. an object will 
form. The task to find how many objects will form in mass 
range (M, M + dM) is converted to the problem of tracing 
the fraction of trajectories passing through the barrier. 

A quantity used to represent the smoothing scale in 
stead of the mass M is the variance of the smoothed field 



S(M) = o 2 {M) = {\S M \ 2 } = Y.( 
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exactly a Brownian random walk. If we define Q(8, S) as the 
number density of trajectories at S within (8,5 + dS), the 
Brownian random walk satisfies a simple diffusion equation 
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and S = 0, S(S) = 0. In absence of barrier, we have solution 
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According to IChandrasekharl il943l) . a trajectory S(S) 
reaches the barrier <5 C at S has equal probability to walk 
above or below the barrier, therefore the solution of Eq. [3] 
with an absorbing barrier boundary condition is 



Q(S,S,8 C 
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The probability of a trajectory absorbed by the barrier 8 C 
must equal to the reduction of trajectories survived below 
the barrier in interval (S, S + dS), 



fs(S,5 c ) = --ig 



QdS 



Substituting Eq.|3] and |S] into the above equation gives 
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which is the fraction of mass associated with halos in the 
range of S and consequently M. So the comoving number 
density of halos of mass at epoch z is simply 
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This is the well-known Press-Schechter mass function. 



(9) 



where St and Wm (k) are the Fourier transform of S and the 
window function Wa/(V) respectively. The smoothed density 
fluctuation can be written as 



5(S) = S(M) =^2S k W M (k) , 



(2) 



which actually tells us that 5(S) is the sum of 8k weighted 
by the window function WmQz). If the smoothing scale R is 
sufficiently large, S and 8(S) will be zero. Once we decrease 
the smoothing scale R, since the window W is a function of 
R, the weighting to Fourier modes of 8 will change. Natu- 
rally the feature of 8(S) trajectories depends on the weight- 
ing pattern of Fourier m odes, i.e. propertie s of the window 
function (see examples in lBond et aljfl99ll) . 

If the window function is sharp in k-space (a top-hat 
function spanning from k — to k ~ 1/R), the increment 
8(S+dS) — 8(S) of a step from S to S + dS comes from a new 
set of Fourier modes in a thin shell of (k, k + dk). Phases of 
8k are uniformly distributed in [0, 2ir], the sum ~Y^~ dk ^fc is 
a random G aussian variable and uncorrelated with previous 
increments feond et al.lll99lt lLacev fc Coldll993l) . This is 



3 THE FRACTIONAL EXCURSION SET 
THEORY 

3.1 motivation 

It is clear that the validity of the Brownian random mo- 
tion prescription to the trajectory of S(S) is guaranteed by 
the sharp k-space filtering. Lack of correlation between the 
new increment with any previous steps delimits the Markov 
nature of the Brownian motion. In context of structure for- 
mation, it means that the formation of halos at small scales 
is not correlated with the density fluctuation smoothed at 
large scales, henceforth halo formation is completely inde- 
pendent of environment. 

If we choose a different smoothing window function such 
as a Gaussian or a top-hat in real space, S(S + dS) contains 
the same set of 8k as S(S) though in the summation each 
Fourier mode is weighted differently by the window func- 
tion. In this circumstance 8(S) is apparently correlated with 
earlier steps, which can not be described by the Brownian 
random walk formalism any longer. In general there is no 
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analytical solution to these types of walks with correlation 
l|Bond et alJll99ll) . 

Recently with the emergence of high resolution sim- 
ulations, it has been revealed that the formation his- 
tory and properties of halo, especially of small mass, are 
modulated significantly by halos' large scale environment 
JSheth fc TormenlbOoilGao et al.ll2005l lHarker et alj|2006t 
IWechsler et al.ll2006h . Therefore there must be considerable 
influence from the mass accretion at large scales on the am- 
plitude of density fluctuation smoothed at small scales, i.e. 
S(S) is correlated with 5(S' < S) even the window function 
is sharp in k-space. 

Either mathematically or physically, the walk of (5(5) of 
a realistic density field is some kind of random motion with 
correlated steps, which is obviously not a Brownian random 
motion, rather, is partly random and partly deterministic. 
Walks like this with "fractional" randomness, fortunately, 
are objects that the fractional Brownian motion (FBM) is 
designed to score. 



3.2 the fractional Brownian motion 

The FBM is a generalizati on of the normal Brown i an ran - 
dom walk introduced by iMandelbrot fc van Nessl l)l96Sl) . 
FBM, though not well-known in astronomy community, has 
been widely used to model geometry and growth of many 
types of rough surfaces in nature like mountain terrain, 
clouds, percolation and diffusion-limited aggregation. Inter- 
estingly it als o finds its application in financial market (c.f. 
lMea,kinlll99^. 

Formally, with index a (0 < a < 1), a FBM is defined 
as a random process X(t) on some probability space such 
that: 

(i) with probability 1, X(t) is continuous and X(0) = 0; 

(ii) for any t > and h > 0, the increment X(t + h) — X (t) 
follows a normal distribution with mean zero and variance 
h 2a , so that 



P (X(t + h) - X{t) < x) 



h~ 



2 /2h 2 



du . (10) 



If a = 1/2, FBM backs to the normal Brownian motion (c.f. 
lFededll988t) . 

The index a is named the Hurst exponent, which is 
used originally in the rescaled range analysis (R/S analysis) 
to portray scaling behaviors of time series. It has strong con- 
nection with the fractal dimensions of time series o r spacial 
struc tures, but the exact relation is case dependent dMeakinl 
1998). Here, a tells us how strongly correlated the step in- 
crement is with previous steps. The trajectory of a FBM 
with smaller index a is more noisy than that of a FBM with 
higher index, so sometimes a is called the roughness expo- 
nent. 

It is very interesting that the FBM has infinitely long- 
run correlations. For instance, the past increments X(0) — 
X(—t) are correlated with future increments X(t) — X(0): 
as X(0) — 0, the correlation function of the "past" and 
"future" is 



_ (-X(-t)X(t)) _ 2a _! 

c{t) ~ 2 1 



(ii) 



which is invariant with the "time" t and only vanishes when 



a = 1/2! This is an impressive feature of FBM, which leads 
us to classify FBM into two types: 

(i) persistence FBM with a > 1/2, which means that an 
increasing trend in past will result in an increasing trend in 
future for arbitrary large t, i.e. a positive feedback process; 

(ii) anti-persistence FBM with a < 1/2, which refers to 
an increasing trend in past will lead to a decreasing trend 
in future, i.e. a negative feedback. 

It might help understanding characteristics of FBM to 
know the generation methods of FBM. To simulate a 1- 
dimensional FBM, the simplest method is 

*" 1/2 G( S ) 



X(t) 



t—1 

s=t — n 



(12) 



r(a + l/2) 

in which G(t) and G(s) are uncorrelated random numbers 
extracted from a normal distribution with zero mean and 
unity variance, n is a practical cut-off number which shall 
be as large as possible. 

To generate a (d + l)-dimensional surface of FBM by 
Fourier transformation, in the first instance we place a grid 
in Fourier space, and fill the grid with complex numbers S(k) 
with Gaussian distributed amplitudes and random phases. 
Spatial correlation is introduced by 



S'(k) = k- (a+d/2) 5(k) 



(13) 



Then Fourier transformation of the random field <5'(k) will 
give a self-affine surface modelled by FBM. 

3.3 the fractional excursion set theory 

The number density of trajectories Q a (8,S) of fractional 
Brow nian moti on with index a obeys the diffusion equation 
fc.f. lLutzll200Hl . 



d 2 Q 

a 



dS ~~ dS 2 ' 
which has solution in absence of barrier 
S~ a ( 8 2 
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Apparently the d istribution of 5(S) at S is still a Gaus- 
sian, the argument of IChandrasekhail (1 1943ft shall be valid, 
thus the solution under boundary condition of a fixed ab- 
sorbing barrier S c is 



) a (S, 5, S c 



-S 2 /2S 2 
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After a straightforward and tedious calculation, the halo 
mass function is 

dlntr 



dlnM 



dU 



with the kernel 



4q 8c 

2^cr 2 



■ exp 



2a 4a 



(17) 
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It is easy see that it is the PS function when a — 1/2. For 
comparison, we reproduce the kernel of the Sheth-Tormen 
mass function here, 
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Figure 1. Ratios of the new mass function of different index a 
to the PS formula. The ratio of the ST function to PS function is 
plotted in thick solid lines, which is in good agreement with the 
new mass function of a 0.435 in high mass regime. The arrow 
indicates the region where all curves cross the line ///ps = 1- 

where A = 0.3222, a = 0.707 and p = 0.3. 

f a of different Hurst exponent a in comparison with PS 
and ST formulas are displayed in Fig. 1. The mass function 
from persistence FBM is very different with that of anti- 
persistence FBM. It appears that the ST function is in good 
agreement with our new mass function of index a = 0.435 
in large mass regime of lna -1 >~ 0.3 beyond which the 
mass function is very sensitive to the choice of a. Consider- 
ing the fact that ST mass function is obtained from fitting 
to simulations and has good accuracy in large mass regime, 
an immediate conclusion is that trajectories S(S) in our uni- 
verse is actually anti-persistent. 

In small mass regime, the dependence on a of f a is 
relatively weak. If a < 0.5 the new mass function predicts 
less number of halos than the PS formula by ~ 10 — 20%, 
but up to ~ 30% more than what ST function gives. It 
is very difficult and unreliable to resolve halos with mass 
lower than 10 8 Mq in present days' simulations, we have to 
leave it to future to tell which mass function is better in 
small mass regime. A quick check indicates that f a has very 
different shape with ST formula at In cr _1 <~ 0.3, we can 
only achieve good fit to fsr in range of -0.5 < In cr" 1 < 0.3 
with a ~ 0.35. 



4 DISCUSSION 

The fractional Brownian motion of index a is introduced 
to construct the fractional excursion set theory. The new 
mass function computed with the theory is analytical and 
simple, of which the PS mass function is only a special case 
of a = 1/2. Comparison with the ST function nurtured by 
N-body simulations demonstrates the excellent performance 
of the new mass function. 

In Fig. 1 it is observed that high mass halo abundance 
is very sensitive to the value of a, the high mass halo abun- 
dance observed can be potentially a very powerful tool to 
detect the non-Gaussianity of the initial density fluctuation 



field: non-Gaussianity will change the correlation between 
walking steps of 5(S) and therefore modify the a effectively. 

The success of applying FBM formalism to model struc- 
ture formation is attributed to the inclusion of the corre- 
lation between density fluctuations at different scales. The 
correlation strength characterized by the Hurst exponent a 
could be resulted from properties of window function and 
the intrinsic correlation of the cosmic density field. We know 
that a non-sharp filtering in k-space will induce correlation 
jBond et alJll99lT) . but are unclear how a changes with fea- 
tures of the window function. More of interests is the re- 
lation between a and the power spectrum of density field. 
The generation method Eq. 1131 provides some clues, how- 
ever there is the complication that the scaling of trajectory 
8(S) is founded relative to the variance a 2 , not the physical 
scale R. Numerical experiments with scale free simulations 
shall be able to improve our understanding effects of window 
function and power spectrum on a. 

In this work only the mass function is computed. In 
principle, the fractional excursion set theory may have many 
applic a tions, for example , th ose works of iLacev fc Cold 
l|l993ft . lMo fc White! lll996h and lZhang fe Huii J2006h can all 
be revisited with FBM. Since a denotes the correlation of <5 
at different scales, the subsequently calculated halo merger 
rate and merger history is marked with the stamp of large 
scale environment on halo formation at small scales. We 
might be able to explain the halo clusterin g dependence on 
halo formation histo ry and environment llGao et alJl2005l : 
IWechsler et"al"ll2006l) . 

The new halo mass function is obtained assuming spher- 
ical collapse. To improve the accuracy of the model, ellip- 
soidal collapse has to be taken into account. The poor per- 
formance of Eq. 1181 of a ~ 0.435 in low mass regime (see 
Fig. 1) is very likely due to our simplification of adopting the 
spherical collapse model. Essentially to calibrate the ellip- 
soidal collapse, we replace the fixed barrier 8 C wi th a moving 
barrier B(S) as in lSheth fc Tormenl (1200211 and lDel PopoTqI 
(2006313), and then solve the diffusion equation with the 
new boundary condition. Technique details and compari- 
son with simulations will be presented elsewhere (Fosalba 
& Pan, in preparation). 
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